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We report ab initio relativistic correlation calculations of potential curves and electric dipole tran- 
sition moments for ten low-lying electronic states, effective electric field on the electron, hyperfine 
constants and radiative lifetimes for the 3 Ai state of cation of the heavy transition metal fluoride 
HfF + , which it is suggested to be used in experiments to search for the electric dipole moment of 
the electron. It is obtained that HfF + has deeply bound X S + ground state; its dissociation energy is 
D e — 6.4 eV. The 3 Ai state is obtained as the relatively long-lived, with lifetime equal to about 0.4 
s., first excited state lying about 0.2 eV higher. The calculated effective electric field, E c g = Wd|fi|, 
acting on an electron in this state is 5.84 x 10 24 Hz/e cm. The obtained hyperfine constants are 
A n = -1239 MHz for the 177 Hf nucleus and A\\ = -58.1 MHz for the 19 F nucleus. 



INTRODUCTION. 

The search for the electric dipole moment (EDM) of 
the electron (d e or eEDM below) remains one of the 
most fundamental problems in physics. Up to now 
only upper limits for \d e \ were obtained. The tight- 
est bound on cL was given in the experiment on the 
atomic Tl beam [l| , which established an upper bound of 
\d e \ < 1-6 x 1CP 27 e cm (e is the charge of the electron). 
Molecular systems provide a way to get much enhanced 
sensitivity, since the effective intramolecular electric field 
acting on electrons in polar molecules can be five or more 
orders of magnitude higher than the strongest field avail- 
able in the laboratory @, H, 0] . 

The new generation of the eEDM experiments, employ- 
ing polar heavy-atom molecules, is expected to reach sen- 
sitivity of 10~ 30 — 10~ 28 e cm/Vday (e.g., see |5j]). Their 
results are expected to dramatically influence all the 
popular extensions of the Standard model, in particular 
supersymmetry, even if bounds on the P, T-odd effects 
compatible with zero are obtained (see [1, 0] and refer- 
ences therein). These studies include the beam experi- 
ments carried out on YbF molecular radicals by Hinds 
and co-workers [j| , and the vapor cell experiment on the 
metastable a(l) state of PbOprepared by the group of 
DcMillc and co-workers (sec [8] and references therein). 
New ways of searching for the eEDM, using trapped cold 
molecular cations, were investigated during the last years 
by Cornell and co-workers. The first candidate was HI + 
[13], but subsequent estimate [l(| and accurate calcula- 
tion [llj have shown that this cation has rather small 
E e g (see below). Then some other candidates were con- 
sidered, HfH+ [12], HfF+ [H] etc., having 3 Ai as the 
ground or probably long-lived excited state. Even experi- 
mental study of spectroscopic properties for these cations 
is a very difficult problem, and there are no such data 
measured up to date though they are already required 
to analyze basic stages of the eEDM experiment [13[ . In 
turn, the modern relativistic computational methods can 
now give reliable answers to almost all the questions of 
interest even for compounds of heavy transition metals 



as is in our case. In the paper 14( a theoretical study of 
the required spectroscopic properties, E c g and hyperfine 
structure, of HfF + as a candidate for the eEDM exper- 
iment was performed. In particular, it was shown that 
3 Ai is not the ground state, its radiative lifetime was esti- 
mated as 0.5 sec. In the present paper the above study is 
disscused in more details and, besides, new obtained data 
are given. In particular, the more accurate data for the 
transition dipole moment 1 £ + — ■> 3 Ai and, consequently, 
more reliable values for the radiative lifetimes of the low- 
est vibrational levels of the 3 Ai are presented. Also the 
obtained data for other transition dipole moments can 
be used to study possible schemes of populating the 3 Ai 
state. With using the Dirac-Hartree-Fock calculations it 
was suggested in [l4| that influence of 4/ electrons of Hf 
is negligible for the excitation energy 1 E + — > 3 Ai. Tak- 
ing into account the importance of this conlusion for the 
accuracy analysis, in the present paper the influence of 
correlation of the 4/ electrons is studied here. 



E eB AND HYPERFINE STRUCTURE. 

One of the most important features of such experi- 
ments is that knowledge of the effective electric field E c g 
seen by an unpaired electron is required for extracting 
d e from the measurements. E e g cannot be obtained in 
an experiment; rather, electronic structure calculations 
are required for its evaluation. It is given by the ex- 
pression E c s = Wd\£l\, where Wd is a parameter of the 
RT-odd molecular Hamiltonian that is presented in Refs. 
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where ^ is the wave function for the 3 Ai state, and f2 = 
(^| J -n\^>) = ±1, J is the total electronic momentum, n 
is the unit vector along the molecular axis directed from 
Hf to F, 
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H d = 2d e 
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E is the inner molecular electric field, and er are the 
Pauli matrices. The E e s value for the 3 Ai state was esti- 
mated in the scalar-relativistic approximation by Meyer 
et al. (sec the note added at the end of [I3|) and the 
method used for calculation of E P ff is close in essence to 
that developed by Titov earlier [18| and applied to two- 
step calculations of the PbF molecule ID, 18| . In paper 
[14j a more reliable value of E e g is calculated using the 
advanced two-step technique developed by Titov and co- 
workers later 0,11 , [21], [24] . The hyperfine constants for 
m Hf and 19 F nuclei (A\\ [Hf] and A\\ [F]) for the 3 Ai state 
of HfF + are also calculated. The hyperfine constants are 
determined by expression [23| 
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where Jhf(f) and MHf(p) are the spin momentum and mag- 
netic moment of 177 Hf( 19 F), a.i are the Dirac matrices for 
the i-th electron, and is its radius- vector in the molec- 
ular coordinate system centered on the Hf or F atom. 
When the experimental value A\\ [Hf] is measured, it pro- 
vides an accuracy check for the calculated E e g value. 
Both A 1 1 [Hf] and An [F] values are useful for identifying 
HfF + by its spectrum from other species in the experi- 
ment. 



METHODS AND CALCULATIONS 
The GRECP method. 

When core electrons of a heavy-atom molecule do not 
play an active role (i.e., their relaxation in the molecule 
is negligible) the effective Hamiltonian with relativistic 
effective core potential (RECP) can be presented in the 
form 



H 



Er 
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Hamiltonian ^ is written only for a valence subspace 
of electrons, which are treated explicitly and denoted by 
indices i v and j v . In practice, this subspace is often ex- 
tended by inclusion of some outer core shells for better ac- 
curacy. In Eq. Q, h Schr is the one-electron Schrodinger 
Hamiltonian 
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where Z ic is the charge of the nucleus decreased by the 
number of inner core electrons. U Ef is an RECP (or 
relativistic pseudopotential (PP)) operator that can be 
written in the separable (e.g., see Ref . |2-4f | and references 
therein) or radially- local (semi-local) [25J approximations 
when the valence pseudospinors are smoothed in heavy- 
atom cores. This smoothing allows one to reduce the 
number of primitive Gaussian basis functions required for 
appropriate description of valence spinors in subsequent 
molecular calculations and also to exclude the small com- 
ponents of the four-component Dirac spinors from the 
RECP calculations, with relativistic effects being taken 
into account by j-dependent effective potentials. Con- 
trary to the four-component wave function used in Dirac- 
Coulomb(-Breit) calculations, the pseudo-wave function 
in the RECP case can be both two- and one-component. 
Besides, the generalized RECP (GRECP) operator^ 



271 ] can be used in Eq. (|4]) that includes the radially- local, 
separable and Huzinaga-type [2^ relativistic PPs as its 
components and some special cases. The GRECP con- 
cept was introduced and developed in a series of papers 
(see Refs. [26, 27, 29|, 3(1 31, 32| and references therein). 
In contrast to other RECP methods, GRECP employs 
the idea of separating the space around a heavy atom 
into three regions: inner core, outer core and valence, 
which are treated differently. It allows one to attain the- 
oretically any desired accuracy, while requiring moderate 
computational efforts since the overall accuracy is limited 
in practice by possibilities of correlation methods. 

Two of the major features of the GRECP version with 
the separable correction described here are generating of 
the effective potential components for the pseudospinors 
which may have nodes, and addition of non-local sep- 
arable terms with projectors on the outer core pseu- 
dospinors to the conventional semi-local RECP operator. 
The problem of division by zero appearing in the cases 
of pseudospinors with nodes is overcome in the GRECP 
method by interpolating the corresponding potentials in 
the vicinity of these nodes [1^, 34 1. It was shown both 
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theoretically and computationally that the interpolation 
errors are small enough. That allows us to generate dif- 
ferent potentials for the cases of outer core and valence 
pseudospinors with the same quantum numbers I and j, 
unlike the conventional RECP approach. In turn, the 
non-local separable terms in the GRECP operator ac- 
count for difference between these potentials, which in 
the outer region is defined by smoothing within the inner 
core as is shown in Ref. 35 , 36|, 37 1 and in many cases this 
difference can not be neglected for "chemical accuracy" 
(about 1 kcal/mol or 350 cm -1 ) of valence energies. The 
more circumstantial description of distinctive features of 
the GRECP as compared to the original RECP schemes 
is given in Refs. [38|, |39( . Some other GRECP versions 
are described and discussed in details in Refs. 2^, 27, 3lll . 

The GRECP operator in the spinor representation[26l 
l40j ] is naturally used in atomic calculations. The spin- 
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orbit representation of this operator which can be found 
m Ref. [HED] is more efficient in practice being applied 
to calculation of molecules. Despite the complexity of 
expression for the GRECP operator, the calculation of 
its one-electron integrals is not notably more expensive 
than that for the case of the conventional radially-local 
RECP operator. 

The GRECP simulating interaction of 12 outercore and 
valence electrons of Hf with the explicitly excluded Is 
to 4/ electrons (60 inner core electrons) is used in 20- 
electron calculations of HfF + . 



Freezing the innermost shells from the outer core 
space. 

The "freezing" of innermost shells from the outer core 
space of electrons within the "small core" GRECPs is 
sometimes required because the accuracy of the GRECPs 
generated directly for a given number of explicitly treated 
electrons cannot always correspond to the accuracy of the 
conventional "frozen core" approximation with the same 
space of explicitly treated electrons (without accounting 
for the frozen states). That space is usually chosen as a 
minimal one required for attaining a given accuracy. It 
was noted in Refs. [3, S3| that using essentially differ- 
ent smoothing radii for spinors with different Ij is not 
expedient since the (G)RECPs errors are mainly accu- 
mulated by the outermost from them. In turn, explicit 
treatment of all of the outer core shells of an atom with 
the same principal quantum number is not usually rea- 
sonable in molecular (G)RECP calculations because of 
essential increase in computational efforts without seri- 
ous improvement of accuracy. A natural way out is to 
"freeze" the innermost of them before performing molec- 
ular calculation but this can not be done directly if the 
spin-orbit molecular basis set is used whereas the core 
shells should be better frozen as spinors. 

In order to exclude ("freeze") explicitly those inner- 
most shells (denoted by indices / below) from molecular 
(G)RECP calculation without changing the radial node 
structure of other (outermore core and valence) shells 
in the core region, the energy level shift technique can 
be applied to overcome the above contradiction [26j, l42| . 
Following Huzinaga et al. [281] . one should add the ef- 
fective core operator U{^ z containing the Hartree-Fock 
field operators, the Coulomb (J) and spin-dependent ex- 
change (K) terms, over these core spinors together with 
the level shift terms to the one-electron part of the Hamil- 
tonian (g]): 



Ug f uz = (J-K) [/„,«] + J2 \fn f li)(fn f lj\ (6) 
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where \ f nf ij)(fn f ij\ are the projectors on the core spinors 
fnfij and Entlj ar e their one-electron energies. The B nf ij 
parameters are presented as M\s nf ij\ in our codes and 
M > 1 to prevent occupying the corresponding states in 
calculations; it is usually selected as M > 1 in our stud- 
ies. Such nonlocal terms are needed in order to prevent 
collapse of the valence electrons to the frozen core states. 
They introduce some "soft orthogonality constraint" be- 
tween the "frozen" and other explicintly treated outer- 
more core and valence electronic states. 

All the terms with the frozen core spinors (the level 
shift operator and exchange interactions) can be trans- 
formed to the spin-orbit representation in addition to the 
spin-independent Coulomb term, using the identities for 
the Pij pro iectors (43]: 



l,j=l±l/2 
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(i.e. ^nfij * fij ~t~ B n y/j) j 



P/ = \lmi)(lmi\ . 

mj=— I 

where 1 and s are operators of the orbital and spin 
momenta, \ljmj)(ljmj\ is the projector on the two- 
component spin-angular function Xljmt 1 \lmi){lmi\ is the 
projector on the spherical function Y\ mi - 

More importantly, these outer core pseudotumors can 
be frozen in calculations with the spin- orbit basis sets and 
they can already be frozen at the stage of calculation of 
the one-electron matrix elements of the Hamiltonian, as 
implemented in the MOLGEP code 44|. Thus, any inte- 
grals with indices of the frozen spinors are completely ex- 
cluded after the integral calculation step. The multiplier 
M=30 was chosen in the present molecular calculations 
to prevent mixing the shifted core states to the wave- 
function due to correlations but not to get poor reference 
wavefunction in the initial spin-averaged calculations at 
the same time (as would be for M — > 00). 

In fact, the combined GRECP version, with separa- 
ble and Huzinaga-type terms, is a new pseudopotential 
treating some minimal number of electrons explicitly but 
which already provide the accuracy approaching to that 
of the frozen core approximation. The efficiency of using 
the "freezing" procedure within the GRECP method was 
first studied in calculations of Tl and T1H 0. 

In ten-electron calculations, which are substantially 
less time consuming, 5s and bp spinors of hafnium and 
Is orbital of fluorine are frozen from the states averaged 
over the nonrelativistic configurations 5o? 2 6s 6 6p ' 4 for 
Hf + and 2s 2 2p 5 for F, and not treated explicitly. 
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The SODCI method. 

The spin-orbit direct configuration interaction 
(SODCI) method is well described in papers 45, Hij ]. 
In the current version of the method, calculations 
are carried out in the AS basis set of many-electron 
spin-adapted (and space symmetry-adapted) functions 
(SAFs). The different AS sets of SAFs are collected 
together in accord to the relativistic double-group 
symmetry requirements for the final configuration 
interaction (CI) calculation. In the present study of the 
molecule having the 2S +1 A' leading term, configurations 
from all the symmetry allowed AS sets with S < 2 are 
included into calculations. 

All the possible singly and doubly excited configura- 
tions with respect to some reference configurations are 
generated (see below for choice of the reference config- 
urations). A generated configuration is included in the 
final CI space if its addition to the reference set leads to 
lowering in the total energy by the value more than some 
threshold T. This lowering is estimated with the help of 
the Ak version of perturbation theory [47$ . in which the 
correlation and spin-orbit interaction are considered as 
perturbations and the wavefunction (obtained from the 
CI calculation in the space of the reference configurations 
for all the AS irreducible representations) is taken as a 
zero approximation (see [22L |42| for details). The lower- 
ings in the total energies for the unselected configurations 
are employed for the T=0 threshold extrapolation. The 
generalized multireference analogue [48| of the Davidson 
correction [i^] (full CI correction) is also calculated. 

In the present calculations, those configurations are 
chosen as the reference (main) configurations which give 
the largest contribution to the wavefunction (i.e. have the 
largest square of the absolute value of the CI coefficient 
(Cj), and thus, that C 2 e/ = S/ercf \ C i I 2 = C > where 
C = 0.973 for ten-electron and C = 0.942 for 20-electron 
calculations. The configurations are obtained from re- 
sults of the preliminary SODCI calculations in the rela- 
tivistic double-group symmetry with the large threshold 
T. New SODCI calculation is then carried out with the 
smaller thresholds. The C is taken the same for each 
point on the potential curve [Hoj ] . Davidson and other cor- 
rections estimating contributions for higher than double 
excitations have an essential dependence on C 2 e y. The 
above selection criterion allows one to stabilize these cor- 
rections for different internuclear distances. This is im- 
portant because the reliability of those corrections has a 
significant dependence on these values. 



Basis sets and property calculations. 



listed in the MOLCAS 4.1 library 53[ was used for flu- 
orine. The molecular orbitals are obtained by the com- 
plete active space self-consistent field (CASSCF) method 
[UlEl with the spin-averaged part of the GRECP 
i.e. only scalar-relativistic effects are taken into account 
at this stage. In the CASSCF method, orbitals are sub- 
divided into three groups: inactive, active, and virtual. 
Inactive orbitals are doubly occupied in all the config- 
urations; all possible occupations are allowed for active 
orbitals, whereas virtual orbitals are not occupied. So 
the wave function is constructed as a full configuration- 
interaction expansion in the space of active orbitals, and 
both active and inactive orbitals are optimized for subse- 
quent correlation calculations of HfF + . According to the 
C 2v point group classification scheme used in our codes, 
five orbitals in Ai, four in Bi and B2, and two in A2 irre- 
ducible representations (irreps) are included in the active 
space. In ten-electron calculations, one orbital in the Ai 
irrep (which is mainly the 2s orbital of F) belongs to the 
inactive space. In 20-electron CASSCF calculations, the 
5s and bp orbitals of Hf and Is orbital of F are added to 
the space of inactive orbitals. 

The ten lowest states with the leading config- 
urations 0^+), [...]<%°\P ( 3 Ai,2,3; X A), 

\..\a\a\n 1 ( 3 n - .o+.i,2! 1 n) were calculated. Here the 
<ti orbital is mainly formed by the 2p z orbital of F with 
admixture of the 6p z and 6s orbitals of Hf, a 2 is mainly 
the 6s orbital of Hf with admixture of the 6p z orbital of 
Hf, 5 and tt are mainly the bd orbitals of Hf. 

To obtain the spectroscopic parameters, six points 
listed in table Q] and a point at 100 a.u. on the HfF + po- 
tential curves were calculated for ten lowest-lying states 
in ten-electron calculations and for four states in 20- 
electron ones. The 20-electron calculation is substan- 
tially more time consuming; therefore the remaining six 
states were calculated for only one point, 3.4 a.u., in the 
present study. Comparing the latter calculations with 
corresponding ten-electron ones, the core (5s 2 and bp 6 
shells of Hf and Is 2 shell of F) relaxation and correlation 
corrections to the T e values, called "core corrections" be- 
low and in Table ITTT| were estimated. Only the electronic 
1 S + state is below the 3 Ai one. The radiative lifetime of 
the lowest vibrational levels of the latter with respect to 
its decay to the vibrational levels of the 1 S + state can be 
written (here we are neglecting the rotational structure 
of the considered states) as 



'el( 3 A!,w2) - 3c 3 
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The generalized correlation atomic basis set 
(12sl6pl6dl0/105)/[6s5p5d3/l5] is constructed for Hf. 
The ANO-L (Us9p4d3f)/[4s3p2dlf] atomic basis set 



(8) 

where AE — (-E(3Ai,u2) — -E(i£+.i,i)) is the difference in 
the energies of the electronic- vibrational states ^3^ i v2 
and ^15]+. di ; d is the dipole moment operator. In the 
adiabatic approximation 
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I < ^At^aldl^iu+.rf > | 2 = 
I < X 3 & 1 ,v2{R)\D{R)iY,+ ^i& 1 \Xi i: + vl {R) > |' 

where t>2(-R) is the vibrational wave function of the 
given electronic state and D(R)i^+_,3/^ 1 is the electronic 
transition dipole moment as a function of the internu- 
clear distance. The excited vibrational levels of the 3 Ai 
state can also decay to the lower vibrational levels of the 
same, 3 Ai, electronic state. The radiative lifetime r v ibr 
of this process is determined by the equations similar to 
© and © with replacing 1 S+ by 3 Ai. £>(i?)3 Al ^3 Al 
here is the total molecule-frame electric dipole moments 
for 3 Ai state of the cation, calculated with respect to the 
center of mass. In Tabic Uwc tabulated D 10e (R) for six 
points using the electronic wave functions obtained in the 
ten-electron calculations. In Table HT1 we listed D 20e (R) 
for point 3.4 a.u. using the electronic wave functions ob- 
tained in the twenty-electron calculations. To calculate 
the radiative lifetimes in according to eqs. ([8|) and (J9j) we 
first use the D 10e (R) values and then multiply the ob- 
tained r e i and T v ibr on the corresponding correction fac- 
tors (£> 10e (3.4)/D 20e (3.4)) 2 . The vibrational wave func- 
tions A3 Al V 2{R)i vi{R) ancl electronic-vibrational 
energies were evaluated on the basis of the 20-elcctron 
SODCI calculations. 

The "atomic core" properties A\\ and E c g were calcu- 
lated only for the 3 Ai state at the point 3.4 a.u., which is 
close to the equilibrium distance (see Table Hn]) . Before 
calculating the core properties the shapes of the four- 
component molecular spinors are restored in the inner 
core region after the two-component GRECP calculation 
of the molecule. For this purpose the nonvariational one- 
center restoration method [l8, 3, 21, 2^, 26|, 5|| is ap- 
plied. The main contribution to the electric field, E, in 
eq. [2] gives the nuclear charge of Hf. In the present paper 
the Hf nucleus is modeled by a uniform charge distribu- 
tion within a sphere with radius i? n uci = 6.8fm. Taking 
into accounting that the electric field from electrons leads 
to a small correction (about one per cent (56J) and the 
dominant contribution is going from the inner core elec- 
trons, we consider here only direct contributions of elec- 
trostatic interaction from the closed shells with n = 1 — 4 
with the EDMs of other (valence and outer-core) elec- 
trons. 



RESULTS AND DISCUSSION. 

The calculated potential energy curves of HfF + arc 
shown in Fig. [TJ The potential energy curves are taken 
from ten-electron calculations but shifted (on the energy 
axis) to fit the D e value obtained in the twenty-electron 
calculations and T e values obtained with twenty-electron 
correction at point R = 3. 4a. it. The results of calcula- 
tions for HfF + spectroscopic parameters are presented in 



Table IIIII The first point to note is that the cation is 
deeply bound. The second important result is that the 
1 E + state appears to be the ground one and the 3 Ai state 
is the first excited one. In addition, the excitation energy 
from 1 £ + to 3 Ai is increased from 866 to 1633 cm" 1 af- 
ter including the 5s and bp shells of Hf and Is shell of 
F into the relativistic correlation calculation. Excitation 
energies from 1 S + to other calculated low-lying states 
are also increased. Note also that the values obtained 
in ten-electron calculations for the lowest four states are 
in a good agreement with the purely 20-electron calcu- 
lations when just the core correction described above is 
taken into account. Accounting for correlation and relax- 
ation of the 4/ 14 — shell would also be desirable for better 
accuracy but it is too time consuming computationally. 
Though the 5s and bp electrons are more strongly bound 
than the 4/ electrons (5s and bp orbital energies about 2 
and 1.2 times higher, respectively) one should rather ex- 
pect that correlation and relaxation contibutions of the 
4/ electrons will lead to a substantially smaller change 
of the spectroscopic properties of HfF + . There are sev- 
eral reasons for this: unlike the 5s, bp shells, the 4/ shell 
is localized in essentially different space region than the 
bd and, particularly, 6s, 6p ones (the average radii of the 
5s and bp shells are about twice larger than that of 4/); 
the angular-type 4/— 4/ and 4/— 6s, 6p electron corre- 
lations are suppressed because there are no "relatively 
low-lying" g-states, whereas the lower-lying 4d-states are 
completely occupied. In Table ITVl the results of the Dirac- 
Fock calculations for the transition bd 1 6s 1 — > 5d°6s 2 with 
[Ae]4/ 14 5s 2 5p 6 and [Ae]4/ 14 frozen core are compared 
to the all-electron results. Two series of Dirac-Fock cal- 
culations are performed; in the first one the cores were 
taken being frozen from the 5d 1 6s 1 state, and in the sec- 
ond one they were extracted from the 5d°6s 2 state. One 
can see from Table HVl that the influence of 4/ relaxation 
is about one order of magnitude smaller than that of bsbp 
relaxation. Similar picture is observed when correlations 
are taken into account. The results of the all-electron 
Fock space relativistic coupled cluster (RCC) calculations 
for the transition from the terms of the 5d 1 6s 1 configu- 
ration to the 5d°6s 2 configuration are presented in Ta- 
ble [Vj Three series of RCC calculations are performed; 
In ten-electron calculation the 5s, bp shells are correlated, 
in sixteen-electron calculation the 4/ shells are corre- 
lated, in twenty four-electron calculation the 4/, 5s, bp 
shells are correlated in addition to the 5d,6s,6p ones. So 
difference between the twenty four-electron and sixteen- 
electron calculations give us the correlation contribution 
of the 5s, bp shells, difference between the twenty four- 
electron and ten-electron calculations give us the contri- 
bution of the correlation of the 4/ shell. One can see from 
Table |W] that the influence of 4/ correlation is about one 
order of magnitude smaller than that of bsbp correlation. 
In the framework of the "atom-in-a- molecule" model, the 
5d 1 6s 1 — > bd°6s 2 excitation of the Hf 2 " 1 " fragment gives 
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the leading contribution to the 1 £ + — > 3 Ai transition of 
HfF + , and therefore we do not expect changes more than 
100 cm -1 in the energy for the above transition when 4/ 
correlation and relaxation are taken into account. So we 
do not expect that the order of levels will be changed. 
The data of transition dipole moments from Table U and 
Table |TT] can be used to study schemes of populating the 
3 Ai state. 

The SO splittings of the 3 A and 3 IT states are mainly 
due to the SO splitting of the 5d shell of Hf: 



H? 



(l5d'S5d) 



(10) 



The atomic Dirac-Fock calculation of Hf + gives 
£ 5/2— £3/2 = 3173 cm -1 => a = 1269 cm -1 , where £5/2 
and £3/2 are orbital energies of the 5d$/2 and 5c?3/2 states 
of Hf + . The SO interaction (flU)) averaged over the 3 A or 
3 II states is reduced to 



H 



LS 



A-(L-S), 



(11) 



where A = 1269/2 = 635 cm -1 , and L and S are the 
orbital momentum and spin moment of HfF + . The SO 
interaction (fTTj) leads to splitting between the compo- 
nents of the 3 A and 3 IT states on 1269 and 635 cm^ 1 , 
respectively. It is in good agreement with the splitting of 
the 3 A state calculated in Table ITTll but with the 3 II split- 
ting because of the off-diagonal SO interaction (which is 
not described by expression (JTTJ)) with closely lying 1 A 
and 1 II states. 

The calculated A\\ and E e g for the 3 Ai state are pre- 
sented in Table PVTl In contrast to A|i[F], AiipEf] and 
i? ff[Hf] are not seriously changed when the outer core 
electrons are included into the calculation. This behavior 
of A\\ [F] is explained by the fact that the shells of fluorine 
in the molecule can be considered with good accuracy as a 
closed-shell subsystem (thus having A\\ [F] w 0) in the 3 A 
state, i.e. there is a large compensation of contributions 
from orbitals with different projections of total electronic 
momentum to Am[F]. Therefore, even a small perturba- 
tion (influencing on the shells of F) can seriously change 
the A\\ [F] value. The calculated E c g is large and com- 
parable with the corresponding value for the a(l) state 
of the PbO molecule [22|. Our E e g value has opposite 
sign than that obtained by Meyer et al. 17] in scalar- 



relativistic calculations and 1.34 times larger by absolute 
value. In Table IVIII we are tabulating lifetimes r e i and 
Tvibr for the lowest vibrational levels of the 3 Ai It is dif- 
ficult to obtain this value accurately because of small 
absolute values of both transition energies and, partic- 
ularly, of the transition dipole moments between those 
states, whereas the absolute errors are similar to those for 
other transitions. Moreover, the transition energy and 
dipole moment calculated for the considered transition 
with rather large relative errors, are presented in equa- 
tion (JSJ) as third and second powers, respectively, thus 



TABLE I: Ten-electron SODCI calculations of transition 
dipole moments and molecule-frame electric dipole moments 
for 3 Ai state calculated with respect to the center of mass. 
Axis z is directed from Hf to F. All values in a.u. 



R 



3.1 3.3 3.4 3.5 3.7 4.0 



D(R)s Al ^3 Al 

D(Rh Al ^ 3ni ■ 
D{R) 3 A 1 -, 3 n 0+ 



ur 

10 1 
10 1 
10 2 
10 2 
■10 1 



-1.15 
1.60 
3.12 
7.55 
3.78 
3.34 
3.86 



-1.40 
1.59 
2.94 
6.52 
4.00 
3.90 
3.30 



-1.52 
1.60 
2.84 
6.02 
4.13 
3.87 
3.01 



-1.64 
1.56 
2.77 
5.59 
4.25 
3.93 
2.80 



-1.87 
1.50 
2.69 
4.93 
4.51 
4.09 
2.45 



-2.20 
1.14 
2.70 
4.11 
4.99 
4.50 
2.13 



seriously increasing the relative error for the lifetimes of 
3 Ai. The correction factors (L> 10e (3.4)/L> 20e (3.4)) 2 cal- 
culated from data of Table Q] and Table |TT] for r e i and 
Tvibr are 0.673 and 1.028, respectively. Finally, it should 
be noted that the calculated energy for the 1 £ + — > 3 IIi 
transition is in a good agreement with the pilot experi- 
mental datum by Cornell group 57 1. 
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TABLE II: Transition and molecule-frame electric dipole moments (with respect to the center of mass for R=3.4 a.u.) obtained 
in twenty-electron calculations. Axis z is directed from Hf to F. All values in a.u. 





x s+ 


3 Ax 


3 A 2 


3 A 3 


X A 2 


3 n 0+ 


3 n! 




3 n 2 




'ill 






-1.20 


1.95-10" 2 











7.77-10" 2 


3.55-10" 


■ l 







5.33-10" 


-l 






-1.50 


1.35-10~ 2 





8.74-10" 3 


3.30-10" 1 


2.56-10" 


■2 


8.42-10" 


■3 


8.23-10" 


-4 


3 A 2 






-1.48 


1.93-10" 


2 3.17-10" 2 





2.31-10" 


■1 


-2.10-10 


-1 


2.15-10" 


-1 


3 A 3 








-1.49 


4.69-10" 2 










0.295 









X A 2 










-1.30 





1.73-10" 


■1 


4.17-10" 


-2 


2.14-10" 


-2 


3 n 0+ 












-9.48-10 -1 


8.42-10" 


3 







1.82-10" 


-3 


3 nr 














-1.05 




2.95-10" 


-2 


1.35-10" 


-1 


3 n 2 


















-1.02 




-5.79-10" 


-2 



x ni -i.24 
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TABLE III: Calculated spectroscopic parameters for HfF + 



State R e A 


T e cm' 1 


T e with core 


w e cm 1 


D e cm 1 








correction" 






10-electron calculation 


1 E+ 


1.784 








751 


51107 




1.810 


866 


1599 


718 




3 A 2 


1.809 


1821 


2807 


719 




3 A 3 


1.807 


3201 


4324 


721 




X A 2 


1.814 


9246 


11519 


696 




X- 


1.856 


9466 


11910 


689 




3 n 0+ 


1.854 


9753 


12196 


699 




3 n! 


1.860 


10190 


12686 


687 




3 n 2 


1.856 


11898 


14438 


703 






1.870 


12642 


14784 


679 




20-electron calculation 




1.781 







790 


51685 


3 A 1 


1.806 


1633 




746 




3 A 2 


1.805 


2828 




748 




3 A 3 


1.804 


4273 




749 





See section "Methods and calculations" for details. 



TABLE IV: Transition energies (in cm x ) for the 5d 1 6s 1 - 
5d°6s 2 obtained in Dirac-Fock calculations of the Hf 2+ ion. 



Frozen 


[Ae]4/ 14 5s"5p b 




all electron 


core 






included 


Frozen from 


13 322 


12 621 


12 439 


5d 1 6s 1 








Frozen from 


11 499 


12 264 


12 439 


5d°6s 2 









[10] B. Ravaine, S. G. Porsev, and A. Derevianko, Phys. Rev. 
Lett. 94, 013001 (2005). 

[11] T. A. Isaev, N. S. Mosyagin, A. N. Petrov, and A. V. 
Titov, Phys. Rev. Lett. 95, 163004 (2005). 

[12] L. Sinclair, J. Bohn, A. Leanhardt, P. Maletinsky, 
E. M. R. Stutz, and E. Cornell, in 36th APS Meeting 
Abstracts (APS, Lincoln, Nebraska, 2005), p. 6119. 

[13] E. Cornell and A. Leanhardt (2006), private communica- 
tion. 

[14] A. N. Petrov, N. S. Mosyagin, T. A. Isaev, and A. V. 

Titov, Phys. Rev. A 75, 030501(R) (2007). 
[15] M. G. Kozlov, V. I. Fomichev, Yu. Yu. Dmitriev, L. N. 

Labzovsky, and A. V. Titov, J. Phys. B 20, 4939 (1987). 
[16] M. Kozlov and L. Labzowsky, J. Phys. B 28, 1933 (1995). 
[17] E. R. Meyer, J. L. Bohn, and M. P. Deskevich, Phys. Rev. 

A 73, 062108 (2006). 
[18] A. V. Titov, Ph.D. thesis, St .-Petersburg (Leningrad) 

State University, Russia (1985). 
[19] A. V. Titov, Int. J. Quantum Chem. 57, 453 (1996). 
[20] N. S. Mosyagin, M. G. Kozlov, and A. V. Titov, J. Phys. 

B 31, L763 (1998). 
[21] A. N. Petrov, N. S. Mosyagin, T. A. Isaev, A. V. Titov, 

V. F. Ezhov, E. Eliav, and U. Kaldor, Phys. Rev. Lett. 



[22; 
[23' 

[24; 
[25' 
[26 
[27 
[28; 
[29; 
[30 



[31 
[32 



[33 
[34 
[35 

[36 

[37 
[38 

[39 

[40 

[41 

[42; 

[43 

[44 

[45 

[46' 



88, 073001 (2002). 

A. N. Petrov, A. V. Titov, T. A. Isaev, N. S. Mosyagin, 

and D. P. DeMille, Phys. Rev. A 72, 022505 (2005). 

. Yu. Yu. Dmitriev, Yu. G. Khait, M. G. Kozlov, L. N. 

Labzovsky, A. O. Mitrushenkov, A. V. Shtoff, and A. V. 

Titov, Phys. Lett. A 167, 280 (1992). 

G. Theurich and N. A. Hill, Phys. Rev. B 64, 073106, 1 

(2001). 

W. C. Ermler, R. B. Ross, and P. A. Christiansen, Adv. 
Quantum Chem. 19, 139 (1988). 

A. V. Titov and N. S. Mosyagin, Int. J. Quantum Chem. 
71, 359 (1999). 

A. V. Titov and N. S. Mosyagin, Russ. J. Phys. Chem. 
74, Suppl.2, S376 (2000), [arXiv: physics/0008160]. 
V. Bonifacic and S. Huzinaga, J. Chem. Phys. 60, 2779 
(1974). 

A. N. Petrov, N. S. Mosyagin, A. V. Titov, and 1. 1. Tupit- 
syn, J. Phys. B 37, 4621 (2004). 

N. S. Mosyagin, A. N. Petrov, A. V. Titov, and I. I. 
Tupitsyn, in Recent Advances in the Theory of Chemical 
and Physical Systems, edited by J.-P. Julien, J. Maruani, 
D. Mayou, S. Wilson, and G. Delgado-Barrio (Springer, 
Dordrecht, The Netherlands, 2006), vol. B 15 of Progr. 
Theor. Chem. Phys., pp. 229-251. 

N. S. Mosyagin and A. V. Titov, J. Chem. Phys. 122, 
234106 (2005). 

A. V. Titov, N. S. Mosyagin, A. N. Petrov, T. A. Isaev, 
and D. P. DeMille, in Recent Advances in the Theory of 
Chemical and Physical Systems, edited by J.-P. Julien, 
J. Maruani, D. Mayou, S. Wilson, and G. Delgado-Barrio 
(Springer, Dordrecht, The Netherlands, 2006), vol. B 15 
of Progr. Theor. Chem. Phys., pp. 253-283. 
A. V. Titov, A. O. Mitrushenkov, and I. I. Tupitsyn, 
Chem. Phys. Lett. 185, 330 (1991). 

N. S. Mosyagin, A. V. Titov, and A. V. Tulub, Phys. Rev. 
A 50, 2239 (1994). 

A. V. Titov and N. S. Mosyagin, in Relativistic Effects 
in Heavy-Element Chemistry (REHE) (Berlin, Germany, 
2003), oral report. 

A. V. Titov, Doctorate Thesis, (Petersburg Nuclear 
Physics Institute, Russian Academy of Sciences, 2002) 
[in Russian]. 

A. V. Titov and N. S. Mosyagin, to be published. 

N. S. Mosyagin and A. V. Titov, arXiv.org/ 

physics/9808006 (1998). 

A. V. Titov and N. S. Mosyagin, arXiv.org/ 
physics/0008239 (2000). 

I. I. Tupitsyn, N. S. Mosyagin, and A. V. Titov, J. Chem. 
Phys. 103, 6548 (1995). 

N. S. Mosyagin, A. V. Titov, and Z. Latajka, Int. J. Quan- 
tum Chem. 63, 1107 (1997). 

A. V. Titov, N. S. Mosyagin, A. B. Alekseyev, and R. J. 
Buenker, Int. J. Quantum Chem. 81, 409 (2001). 
P. Hafner and W. H. E. Schwarz, Chem. Phys. Lett. 65, 
537 (1979). 

A. V. Titov, A. N. Petrov, A. I. Panin, and Y. G. Khait, 
program package for calculation of molecular matrix ele- 
ments with the Generalized RECP. 

R. J. Buenker and S. Krebs, in Recent Advances in Mul- 
tireference Methods, edited by K. Hirao (World Scientific, 
Singapore, 1999), pp. 1-29. 

A. B. Alekseyev, H.-P. Liebermann, and R. J. Buenker, in 
Recent Advances in Relativistic Molecular Theory, edited 
by K. Hirao and Y. Ishikawa (World Scientific, Singapore, 



10 



TABLE V: Transition energies (in cm x ) from the terms of the 5c! 1 6s 1 configuration to the 5d°6s 2 configuration obtained in 
RCC calculations of the Hf 2+ ion. 



Transition 




ten-electron sixteen-electron twenty four-electron 


conribution 


conribution 


from 




calculation 


calculation 


calculation 


from 5s5p 


from 4/ 


6si /2 5d 3 / 2 (J = 


1) 


9233 


10899 


9182 


-1717 


-51 


6s 1/2 5d 3/2 (J = 


2) 


8935 


11353 


8817 


-2536 


-118 


6s 1/2 5d 5/2 (J = 


2) 


6251 


8763 


6036 


-2727 


-215 


6s 1/2 5d 5/2 (J = 


3) 


5168 


6991 


4973 


-2018 


-195 



TABLE VI: Calculated parameters A,|[Hf] and A||[F] (in 
MHz) and E cB (in 10 24 Hz/(e • cm)) for the 3 Ai state of 
177 Hf 19 F + at internuclear distance of 3.4 a.u. 



A,|[Hf] A|,[F] 



10-electron calculation 



-1250 



-33.9 



5.89 



20-electron calculation 



-1239 



-58.1 



5.84 



TABLE VII: Calculated vibrational energy levels (in cm -1 ) of 
the 1 S + and 3 Ai states and radiational lifetimes (in seconds) 
of the 3 Ai state. 



, A 1 



E v 
416 368 
1234 1117 
2038 1878 
2821 2631 
3586 3373 



Tel 

0.389 
0.365 
0.343 
0.316 
0.290 



0.184 
0.088 
0.059 
0.046 
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